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We show that the strong couphng of pseudospin orientation and charge carrier motion in bilayer 
graphene has a drastic effect on transport properties of bahistic p-n-p junctions. Electronic states 
with zero momentum parallel to the barrier are confined under it for one pseudospin orientation, 
whereas states with the opposite pseudospin tunnel through the junction totally uninfluenced by the 
presence of confined states. We demonstrate that the junction acts as a cloak for confined states, 
making them nearly invisible to electrons in the outer regions over a range of incidence angles. 
This behavior is manifested in the two-terminal conductance as transmission resonances with non- 
Lorentzian, singular peak shapes. The response of these phenomena to a weak magnetic field or 
electric-field-induced interlayer gap can serve as an experimental fingerprint of electronic cloaking. 



Charge carriers in graphene behave like relativistic 
particles Some of the most intriguing aspects of 
carrier dynamics in this material arise due to chirality, 
i.e. the strong coupling of pseudospin and orbital degrees 
of freedom [2]. Analogs of relativistic electron effects such 
as Klein tunneling [3^-^], and optical phenomena such as 
negative refraction Q, Fabry- Perot resonances [3-[2| and 
the Goos-Hanschen effect [lO|, ujj , provide a platform for 
understanding transport in graphene nanostructures [12I- 
14j . Chirality was also proposed as a vehicle for coupling 
the orbital motion of carriers to the inner valley degrees 
of freedom of graphene 15|, l^] • 

Here we describe new effects in bilayer graphene (BLG) 
that have no direct analogs in optics or in single-layer 
graphene. These effects, which arise due to the chiral na- 
ture of carriers, have dramatic consequences for transport 
through potential barriers. At normal incidence, chiral- 
ity mismatch leads to complete decoupling of states in 
regions of opposite polarity (see Fig|T^,b). Electrons of 
one chirality are confined within the barrier, despite the 
presence of a continuum of available states outside the 
barrier. Conversely, electrons in the outer region scatter 
and tunnel through the barrier as if no localized states 
were available on their way. Thus the barrier acts as 
a cloak for confined states, rendering them invisible via 
both transmission and reflection. This effect, which is 
unique to BLG, leads to a number of intriguing and po- 
tentially useful properties of BLG nanostructures, such 
as tunable confinement and coupling to individual states. 

For oblique incidence, pseudospin decoupling is imper- 
fect, producing transmission resonances associated with 
confined states. However, these resonances are very nar- 
row at near- normal incidence angles, making the outer 
states 'blind' to the confined states at almost all energies 
(see Figs, [it and [2j). This yields a high-fidelity cloak 
effect for all near- normal incidence angles. The effect is 
sensitive to external magnetic field applied perpendicular 
to BLG plane: the Lorentz force curves particle trajec- 
tories, making the orbits normally incident on the bar- 
rier exit at oblique angles. We show that the narrowing 
of transmission resonances at normal incidence is mani- 
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FIG. 1: Cloaking of transmission resonances in ballistic p-n-p 
junctions due to decoupling of states with orthogonal pseu- 
dospins. For a potential step (a), the pseudospin-reversing 
coupling (wavy line) vanishes at normal incidence. This leads 
to formation of confined states under a barrier (b), seen as 
narrow resonances in transmission. The cloak effect illus- 
trated in (c) for a delta function barrier, having one confined 
state which produces a transmission resonance e = (see 
Eq.© and accompanying text). The resonance is narrowed 
and dimmed at near- normal incidence angles (small q). As 
seen in the inset, showing trace along dashed line, transmis- 
sion at all e / eo is uninfluenced by the confined state (the 
cloak effect). 

fested in sharp non-Lorentzian peaks in the two-terminal 
conductance with square-root singularities at the tips 
(see Figs. [2] and [3]). This singular behavior, along with 
its suppression in magnetic field, can serve as an experi- 
mental signature of the cloak effect. 

It is instructive to compare the electron cloak real- 
ized in BLG with optical cloaks of invisibility 0, fisj , 
which employ refractive metamaterial shells to guide 
plane waves around an object. In contrast, in our ap- 
proach, the probing wave is transmitted directly through 
the region containing the cloaked states, made invisible 
by decoupling of opposite chiralities. Crucially, cloaked 
states can either hold particles or be empty. Since elec- 
tron interactions in BLG are nearly pseudospin-blind[2|, 
filling of confined states does not affect the decoupling 
of opposite chiralities which is responsible for the cloak 
effect. For similar reasons, states with different spin and 
valley polarization can be cloaked simultaneously. 
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FIG. 2: The transmission probability (a) and the conduc- 
tance (b) for parabolic barrier U{x) — —^moj^x^. Shown are 
resonances corresponding to the first two harmonic oscillator 
bound states, e = —^huu and —^huu. At normal incidence, 
^ = 0, these confined states are cloaked by the barrier due to 
chirality mismatch with the continuum states outside. The 
conductance resonances feature the characteristic square-root 
singularity 6G oc — i/|^e|. Units: = GoW/i, where W is 
the sample width and i = ^h/2m(jj. A zoom-in shows the 
first resonance overlaid with 10 transmission profiles, with qi 
changing from 0.18 to as indicated by an arrow. 

We also note that scattering on barriers in BLG was 
studied in Refs.[3, 8], where Fabry-Perot resonances at 
oblique incidence were found. Analytic results for trans- 
mission through a square barrier at normal incidence 
were obtained [3], however the cloak effect and its rela- 
tion to confined states was not elucidated. 

In experiments, several factors may modify the picture 
of chirality-induced cloaking of confined states. First, 
disorder in the barrier region, which breaks conserva- 
tion of the parallel momentum contributes additional 
broadening of the resonances, in particular near q = 0. 
However, recent advances in experimental techniques, 
such as using suspended graphene samples [HI, and de- 
position on hexagonal boron nitride [2o| , help to dramati- 
cally enhance the mobility of BLG samples and minimize 
the influence of disorder. Second, the top and back gates 
used to create p-n-p junctions in graphene devices natu- 
rally create a potential difference between the two layers, 
thus opening a gap in the BLG spectrum 21, |22|. The 



effect of gap opening will be discussed below. 

Electronic transport in BLG in the presence of a po- 
tential barrier U{x) is governed by the Hamiltonian[22| 



H = 



U{x) 

2m V dT 



2m V dec ^ 

U{x) 



(1) 



where the barrier is spatially uniform in the ^/-direction, 
and q is the conserved wavevector component parallel to 
it (m ^ 0.04me is the BLG band mass). The cloak ef- 
fect is complete for normally incident electronic states, 
for which q = 0. The problem can be conveniently an- 



alyzed in the eigenbasis of cr^^, ax\^) = In this 

basis the pseudospinor Schrodinger equation decouples 
into two scalar Schrodinger equations: 



e) i>±{x)=0. 



(2) 



For the special case of a potential step (see Fig(T^), and 
for energies between the values of U{x) far to the left 
and far to the right, /7n < e < /7p, plane wave states on 
one side are coupled to evanescent states on the opposite 
side. Thus, despite the availability of free carrier states 
with energy e on either side of the step, the transmission 
is completely blocked due to chirality mismatch. 

Extending this argument, we show that chirality leads 
to the confinement of electronic states under potential 
barriers in BLG. For a general form of potential bar- 
rier /7(x), as shown in FigHb, the two sign choices in 
Eq.(|2]) lead to two qualitatively distinct types of states 
with zero momentum parallel to the barrier {q = 0). The 
solutions for describe continuum states that live 

outside the barrier region, while the solutions for 
describe bound states confined in the inverted potential 
— U{x). Because these confined states completely decou- 
ple from the continuum states at g = 0, they do not show 
up in the normal incidence transmission (c.f. Ref.[3]). 

Crucially, transmission through the barrier in BLG is 
controlled by different mechanisms at normal incidence 
and at oblique incidence. At normal incidence, since 
pseudospin is conserved, the only available mode of trans- 
mission is direct tunneling via the evanescent wave that 
extends through the barrier. The character of transmis- 
sion changes completely for oblique incidence. Extending 
Eq.(|2]) to nonzero we obtain two coupled equations: 

dx'^ 



V^±(x) = 2q 



d^^{x) 



dx 



(3) 



Due to the ipj^/il)- mixing, the confined states acquire fi- 
nite coupling to the continuum, turning into transmission 
resonances. However, since the width of these resonances 
vanishes as q^ at near- normal incidence, the cloak effect 
persists in a finite range of incidence angles with small q. 

We highlight various aspects of the cloaking behavior 
by analyzing two models, a narrow barrier and a wide 
barrier, modeled by delta function and inverted parabolic 
potentials. In the first case, U{x) = A(5(x), A > 0, at 
zero q we get a confined state iIj-{x) oc e~'^°l^l, with 
= Xm/fP and energy eo = \^m/2TP. Scattering states 
for a delta function can be found exactly, by extending 
the standard approach to account for evanescent states 
which appear in the free particle BLG problem. 

We nondimensionalize Eq.(j2j) via e (fi^/2m)e, A 
(fi^/2m)A, and write the free particle wavefunction as 



'ij)±[x > 0) = a±e 



c±e 



e>q\ (4) 



with k = + g^, k = — cp. For a potential with 
inversion symmetry, \J{x) = /7(— x), there are two types 
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of symmetric solutions of Eq.(j2j): those corresponding 
to even and iIj-{x) odd, and vice versa. Taking 

?/^+(— x) = tlj-{—x) = ^tlj-{x) in each case, and 

performing standard matching for the values and deriva- 
tives at X = 0, we find a relation between in and out 
states of the form 6^'^^ = 2;(i^2)C+'^\ where the label 1 
(2) indicates that ^l^-\-{x) is even (odd). The transmission 
amplitude for a plane wave incident from one side is found 
by taking a suitable superposition of these even/odd par- 
ity states, giving t(e, q) = ^{zi — Z2). In this way we find 



X-^ + irX- 



A_ 



irXk X- 



irX^ 



(5) 



where X±^ = A ± 2/^, X±k = A ± 2z/c, r = /kn. At 
small q^ the first term exhibits a resonance associated 
with the confined state, with the energy eo, whereas the 
second term is nonresonant. As shown in Fig(T]3, the 
cloak effect suppresses the transmission resonance e ^ eo 
at near- normal incidence angles. Away from resonance, 
transmission as a function of e and q closely mimics that 
for the Schrodinger problem, Eq.(|2]). This is illustrated 
in Figdt inset. 

A similar behavior was found for a wide barrier, taken 
to be U{x) = — ^mcj^x^. The transmission probability, 
obtained numerically, is shown in Figl2^. For normal 
incidence, this potential supports a family of chirality- 
induced bound states with evenly-spaced energies = 
—fujoin + 1/2). Accordingly, our simulation yields res- 
onances in transmission and conductance peaked near 
these values (Figs J2|3p . The cloak effect is manifested in 
resonance widths vanishing at normal incidence. 

Since the BLG Hamiltonian supports coexisting 
evanescent and plane wave states, the procedure for 
this solution is notably more complicated than that for 
the non-relativistic Schrodinger equation. In particular, 
Eq.(|3]) must be solved twice for each once with an out- 
going plane wave as the initial condition, and once with 
a decaying evanescent wave. The physical solution is the 
linear combination of these two solutions in which the co- 
efficients of the growing evanescent waves on both sides 
are zero. We evaluate the conductance using the relation 



G{e) = GQW / dq\t{e,q)\\ Go 



(6) 



where W is the width of the sample in the ^-direction, 
and TV = 4 is the spin/ valley degeneracy of BLG. The 
integral in Eq.(j6]) runs over the interval —kp < q < kp, 
where /cf is the Fermi wave vector in the leads. 

Unlike the case of a delta function barrier, where reso- 
nance appeared as a weak feature on top of a large back- 
ground transmission, transmission at zero q is fairly low 
for the parabolic barrier. This is so because in the cloak- 
ing regime transmission is dominated by direct tunnel- 
ing without coupling to confined states. Since tunneling 
is exponentially small for a wide barier, transmission is 
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FIG. 3: The effect of interlayer gap on ballistic conduc- 
tance for a parabolic barrier, U{x) = —^muj'^x^. Cloaking 
manifests itself in abnormally sharp resonances positioned at 
e = —huj{n + |), with n — 0,1,2,3 for a gapless system. 
The sharp cusps observed in the gapless case, Eg. fTT)) . are 
rounded by a nonzero gap. Results for gap values = 0, 
2/kj, and Ahw are shown. Zoom- in on the resonances at 
e = —^huj^—^huj (inset) indicates that an overall suppres- 
sion of conductance, estimated from the max/min ratio, has 
little effect on the fringe contrast. 

dominated by resonances. The reduction of direct tun- 
neling compared to resonant transmission is clearly seen 
in simulation (see FiglJjD inset). The nearly ffat (black) 
line is the background transmission, calculated for g = 0, 
which is due to direct tunneling under the barrier. The 
asymmetric lineshape for nonzero which arise due to 
interference of the direct and resonant transmission path- 
ways, can be described by the Fano model(23|. 



^Fano(e) = 



l + a2 (e-e,)2+^2' 



(7) 



where 7 and a are the energy width and 'asymmetry pa- 
rameter' that depend on the specifics of the system. 

Higher order resonances display more symmetric (ap- 
proximately Lorentzian) profiles which can be explained 
by large values of the parameter a in Eq. ([7]). Large a 
values indeed follow from an estimate based on the gen- 
eral relation derived in Ref.[2^, a = |ro/to|, where tq and 
to are the background refiection and transmission ampli- 
tudes in the absence of the resonance. For a wide barrier, 
the problem can be described by the WKB approach, giv- 
ing transmission which is exponentially small, and thus 
a large a value. In particular, for a parabolic barrier we 
have \ro/to\ = exp {ire/hw). Hence for values of e above a 
few times ficj, the transmission amplitude near resonance 
is well- approximated by the Breit-Wigner model. 



t{e,q) 



l[q)^i{e- e^[q)y 



(8) 



where the parameters e^{q) and ^{q) describe the reso- 
nance energy and width as a function of wave vector q. 
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What is the origin of the pecuhar peak shapes seen in 
conductance? To elucidate their relation with Lorentzian 
peaks in transmission, Eq.(|8]), we examine the quantity 
SG{e) = G{e) - G(e*(0)). Using Eq.©, we write 



6G{e) = GoW 



\t{e,{0),q)\')dq. (9) 



Since the difference of the two terms under the integral 
quickly goes to zero for 7(g) ^ |e — e*(0)|, near the res- 
onance the integral is dominated by small q. 

Using quadratic model dependencies e^{q) = e*(0) + 
aq'^, 7(g) = /3q'^ which are valid at small we can sim- 
plify the expression for SG as follows 



5G{e) 



GoWP^ 



^^^^^ ^ a 



where Se = e— e*(0), u = \^q. The parameter a provides 
a relative measure of the rates at which the resonance 
moves and widens as a function of q. Integration can be 
done using the partial fraction decomposition, 



{a + i)Se 



+ {5e — au'^)'^ + i{au^ — de) 



ex., 



(10) 



and the identity au'^+h ~ 7r(a6)~-^/^, where for com- 
plex a and h the branch of the square root corresponding 
to Re \Jhl a > should be used. This gives 



6G{e) = -GoW- 



30 
2 



f sgn((5e) 




(11) 



(j) = arctan(Q;). This analysis predicts a square root sin- 
gularity near the tips, SG{e) oc ^|Je|. The values a and 
P are determined by the details of the barrier potential, 
and lead to four possible singular lineshapes [25]. Our 
numerical results confirm this analysis: peaks in trans- 
mission are approximately Lorentzian (see Fig. [2] inset), 
whereas the conductance peaks feature sharp cusps ap- 
proximately described by 6G{e) oc Y^|(5e|, see Figs. l2|3[ 

As discussed above, the cloaking of confined states be- 
comes imperfect in the presence of magnetic field, since 
the cloaking condition q = is violated by trajectory 
bending. This effect can be incorporated in our anal- 
ysis via a momentum shift q ^ q -\- {eB/ch)x, giving 
6q = {eB/ch)L for a change in momentum of a particle 
moving across the barrier, where L is the barrier width. 
The effect of a weak field can be mimiced by introduc- 
ing a cutoff Sq in the integral ([9j), leading to rounding 
the square-root cusps on the scale Se ~ B^L^. Thus, 
uncloaking of confined states in weak fields will manifest 
itself in suppression and rounding of conductance peaks. 

It is straightforward to incorporate the effect of in- 
terlayer gap opening in the analysis [25]. As shown in 
Figini the opening of a gap weakly modifies the positions 
and shapes of the resonances, but does not destroy them. 



In particular, the fringe contrast, as assessed by the 
max/min ratio, is relatively unaffected by the gap (see 
Figl3] inset). Therefore the resonances in conductance 
provide a robust signature of confined states. The peak 
shapes are singular when these states are cloaked, and 
smeared when they are uncloaked by applying a magnetic 
field or an inter layer gap. These features provide a char- 
acteristic signature of cloaking and uncloaking, which we 
expect to be readily observable in experiments. 

The phenomena described above originate from decou- 
pling of chiralities for normally incident states leading to 
formation of confined states within the barrier. Because 
such behavior is diametrically opposite to that studied in 
single layer graphene, it might be called 'anti-Klein tun- 
neling'. However, since the term 'Klein tunneling' is in 
a sense a misnomer describing a non- tunneling behavior, 
we prefer to avoid introducing distinctions based on the 
degree of kleinness. 

To summarize, chirality mismatch of states inside and 
outside a ballistic p-n-p junction results in cloaking of 
states confined inside the junction. The cloak effect 
is perfect at normal incidence angle, and close to per- 
fect at near- normal incidence. The confined states man- 
ifest themselves as resonances in the ballistic conduc- 
tance with characteristic non-Lorentzian lineshapes with 
square-root singularities due to cloaking. These singular 
resonances, which are smeared by weak magnetic fields, 
can serve as a hallmark of the cloak effects in transport 
measurements. 

This work was supported by Office of Naval Research 
Grant No. N00014-09-1-0724. 



APPENDIX A: SINGULAR SHAPES OF THE 
CONDUCTANCE PEAKS 

As shown in the main text (see Eq. (pT]) and accompa- 
nying discussion) , the shape of conductance peaks is sen- 
sitive to the relative rates at which the transmission res- 
onances shift and widen as a function of the parallel mo- 
mentum q. The different types of behavior, illustrated in 
Fig. m can be understood as follows. In the limit a <C /3, 
the position of the resonance is approximately stationary 
as a function of q, producing a symmetric lineshape (since 
the transmission profile is symmetric). In this case, SG{e) 
is negative, forming a cusp at e = [Fig. HJ^a)]. As a in- 
creases, the contributions of transmission with different 
q values move to one side of the resonance, resulting in a 
conductance hneshape that is asymmetric about e*. For 
a sufficiently rapidly shifting resonance, at the critical 
value ac = P/V^^ the leading term for SG{e) on one side 
of the resonance vanishes, giving a flattened behavior SG 
vs. e [Fig. m^b)]. Increasing a above ac changes the sign 
of SG{e) [Fig. m^c)]. For a the resonance position 

moves fast as a function of q resulting in nonzero G{e) 
on only one side of the resonance [Fig. HJd)]. 
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ax = ±1, we obtain a generalization of Eq. (|3j): 
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FIG. 4: Schematic for the conductance peak shape G{e) for 
different values of the parameter a describing the dispersion of 
the quasibound state, e^{q) — e^-\-aq^, Eq. (|lip . Small a ^ /S 
corresponds to a cusp with a square-root profile (a). For a 
near the critical value ac = the peak flattens out on one 

side (b). For a > the peak transforms into a monotonic 
lineshape (c), (d). The results for positive a (shown) and for 
negative a are related by symmetry 6e ^ —6e. 
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FIG. 5: Transmission and conductance of the gapless system 
(top), and gapped systems with Eg = huo (middle) and Eg = 
4huj (bottom). Axes labels are the same in Fig. [2] Shown 
are resonances correspond to the harmonic oscillator states 
e = — l/kj, — — —^huj. As the gap Eg increases, it 
starts to dominate the resonant tunneling coupling, producing 
a rounding of the suare-root singularities at the peak tips. 



APPENDIX B: THE INTERLAYER GAP EFFECT 

The generalization of our problem to BLG in the pres- 
ence of non-zero transverse polarization that opens up an 
inter layer gap Eg is described by adding a term ^Egas 
to the Hamiltonian ([1]), i^gap = H i^g^s- Writ- 
ing the Schrodinger equation in pseudospin components 



2 , "^rn,^^, , 



^ d ^ mEg 
dx h? 



(12) 

From the coupled equations Eq. (fT2|) . we see that for 
Eg ^ and at g = 0, there is a coupling between the 
confined and deconfined states ?/^+ and i/j- of order Eg. 
This coupling limits the width of a resonance in transmis- 
sion by ^{q) (X when q = 0^ which sets the minimum 
cusp width for conductance peaks to be of order E'^ (see 
Fig. [5]). We also note that there is a line of transmis- 
sion zeros away from g = 0, which corresponds to the 
situation where the right hand side of Eq.(|3]) is zero. If 
we approximate (^)~^ as A = h/\^2me^ we find a curve 



described by q£ 



Eg/hoj 



which shows that the trans- 



mission zero gets closer to q = for increasing e. Also, 
the curve collapses to q = when Eg goes to zero, as 
expected from the cloak effect. 
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